Quantum computation with abelian anyons on the honeycomb lattice 
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We consider a two-dimensional spin system that exhibits abelian anyonic excitations. Manipu- 
lations of these excitations enable the construction of a quantum computational model. While the 
one-qubit gates are performed dynamically the model offers the advantage of having a two-qubit 
gate that is of topological nature. The transport and braiding of anyons on the lattice can be per- 
formed adiabatically enjoying the robust characteristics of geometrical evolutions. The same control 
procedures can be used when dealing with non-abelian anyons. A possible implementation of the 
manipulations with optical lattices is developed. 



Topological models address in the most sufficient way 
the problem of decoherence, which is the main obstacle 
in the realization of quantum computation P, Q IS ■ 
This is achieved by encoding information in the sta- 
tistical properties of non-abelian anyons. Hence, infor- 
mation remains unaffected by temperature or local per- 
turbation as long as the nature of the anyons does not 
change. Subsequently, information can be processed by 
braiding anyons, an evolution which is topologically pro- 
tected from control errors. Nevertheless, due to their 
exotic nature these models are hard to realize physically. 
To approach this problem we construct a simple model 
for quantum computation with topological characteris- 
tics that can be implemented easier than the fully topo- 
logical quantum computer. The simplest case is to iden- 
tify the topological phase arising from the statistics of 
abelian anyons with a two-qubit gate 5] . In this scenario 
qubit encoding and one-qubit gates can be performed in 
the usual dynamical way. 

The purpose of this letter is twofold. First, we propose 
a physically realizable model for the generation and de- 
tection of abelian anyons. The latter can be subsequently 
used to perform quantum computation that has a topo- 
logically protected two-quit gate. Second, we develop the 
test ground for the generalization to the more demanding 
non-abelian case where topological quantum memory or 
computation can be performed. Indeed, the main ingre- 
dients of the abelian case, such as the creation of anyons, 
their braiding and their subsequent measurement can be 
carried through to the case of non-abelian anyons. 

While the ideas here can be applied to a variety of dif- 
ferent setups 0, we shall employ a particular physical 
model, namely a two dimensional spin-1/2 system where 
the spins are located at the vertices of a honeycomb lat- 
tice |a] . The spins are set to interact with each other via 
the following Hamiltonian 
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where "x-links", "y-links" and "z-links" label the three 
different link directions of the honeycomb lattice as de- 
picted in Figure n It has been shown Q that this model 
exhibits two main distinctive phases. For | Jj;| + |Jj,| > 



l-^^l, l-'^yl + \Jz\ > \Jx\ and |J^| + \Jx\ > \Jy\ a gap- 
less sector appears that supports non-abelian excitations, 
while the violation of any of the three inequalities leads 
to a gap. Here we shall focus on the latter phase for a 
concrete choice of couplings, namely \Jz\ > \ Jx\ + \ Jy\- 
As we shall explicitely see here this phase is equivalent 
to Kitaev's toric code which exhibits three distinguished 
species of particle excitations above the ground state 0, 
namely a fermion, e, and two bosons, the chargeon, e, 
and the magnon, m Interestingly, exchanges between 
different particles reveals their anyonic character. 

To identify the properties of the abelian phase in the 
honeycomb lattice it is instructive to employ perturba- 
tion theory 8]. Nevertheless, the discrete qualitative 
properties of the model will stay the same even beyond 
the perturbative approximation 0. Consider the cou- 
pling regime where JxiJy,Jz > and Jz ^ Jx,Jy The 



dominant z-interaction has as lower eigenstates 



I tt) and [j = III) for each z-link. The Pauli opera- 
tors that act on these states are given by = a'^a^. 
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a^l 1(7^, where 



stands for equality in the space spanned by the fy and 

vectors. Considering up to fourth order in perturba- 
tion theory, the Hamiltonian becomes 
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where an irrelevant constant term has been omitted. The 
summation in Q is over all plaquettes of the lattice, 
while the lower indices of the Pauli operators denote the 
position of the effective spins around the plaquette, p, 
starting from the left and moving clockwise (see Fig- 
ure 

This Hamiltonian is unitarily equivalent to the more 
familiar Kitaev's toric code 1], but here we shall work 
with Hamiltonian jSJl. Before studying the particle ex- 
citations we first determine the ground state, | 0) . Since 
Qp — Qp and Qp = 1, the eigenvalues of Qp are ±1 
and since [Qp,Qp'] — for all p, p' , the lowest energy 
state has to satisfy Qp | 0) = | 0) for all p. Excitations of 
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Hamiltonian (0) are given by relaxing this requirement 
at some plaquettes. For periodic boundary conditions we 
have Yip Qp = 1 at all times, hence, excitations should 
appear in pairs. For open boundary conditions single ex- 
citations are possible, where the other part of the pair 
can be considered to be outside the system. Neverthe- 
less, for manipulations away from the boundary we may 
assume that the generated particles come in pairs. 
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FIG. 1: The honeycomb lattice with the three different types 
of links X, y and z. The distribution of the four spin inter- 
action, Qp, is indicated. Possible positions of the excitation 
pairs produced by rotations of effective spins. 

Consider a rotation ct^ applied on a certain z-link. As 
anticommutes with the Qp interactions of the plaque- 
ttes on its left and on its right it will increase the energy 
of these plaquettes by AE = 2 Jeff. Similarly, if we per- 
form a rotation ct^, then the energy of the plaquettes 
above and below increases by AE. If a rotation ia^ is 
performed, then all four neighboring plaquettes obtain an 
energy increase of AE each. Hence, the low lying energy 
eigenstates are elementary pairs of excitations given by 

\Z)=a]\0) and \Y)=d-]\0). (3) 

These excitations live in appropriate neighboring plaque- 
ttes as indicated in Figure ^ Moreover, one can combine 
such pairs of excitations in a chain of and rota- 
tions to create two excitations at the end points of an 
arbitrary string s. Similarly, one can define a pair of 
fermionic excitations, \X), as 

\X[s]}^l[^a^\0), (4) 

jes 

where s is the path along which tr^ rotations are per- 
formed in such a way as interweaving of the constituent 
CT^ and paths is achieved. The eigenvalues of the en- 
ergy of the Y and Z excitations have an energy gap above 
the ground state given by AE = 4Joff, while for the X 
excitation the energy gap is double. Note that, the eigen- 
value of the energy of each excitation is independent of 
the length of the string, s, as the contributions to the 
energy come only from its two end-points. 

From these definitions one can easily deduce the fu- 
sion rules XxX = YxY = ZxZ^l, XxY = Z, 
Y X Z = X and Z x X = Y, where 1 denotes the vac- 
uum. Indeed, two successive identical Pauli rotations (X, 



Y on Z), that correspond to superposing two pairs of the 
same excitations, will give the identity operation. The 
generation of an X-particle out of the fusion of Y- and 
Z-particles follows straightforwardly from the property 
(jJ'cr^ = ia^. Furthermore, from definitions ||2J) and Q 
one can deduce the statistics of these particles. Exchang- 
ing two X-particles results in a complete overlap of their 
strings at 2{2k + 1) sites, where k is an inteeger, giving 
eventually {ia^Y = — f as an overall factor that reveals 
their fermionic character. To probe the anyonic char- 
acter we circulate an X-particle around a single Y- or 
Z-particle. This loop that corresponds to two consecu- 
tive exchanges gives a string operator that anticommutes 
with the string of the Y- or Z-particles, due to the anti- 
commuting properties of the Pauli operators at the site 
of their intersection. Hence, a minus sign will be pro- 
duced revealing that the X and Y oi Z particles behave 
as anyons with respect to each other with statistical angle 
9 = 7r/2. Note, that due to the undetermined structure 
of the strings for the Y and Z particles it is not straight- 
forward to deduce their statistical properties. 

In order to use this setup to perform quantum com- 
putation we encode the states of a qubit in the space 
of excitations of the system, namely |0l) ~ \X) and 
= \ Y) (or = \Z) depending on the position 
of the excitation). To manipulate these states we place 
them in neighboring plaquettes. A one-qubit rotation, 
e.g. of the I Oi) state, can then be performed by an ap- 
propriate z-rotation, e~*'^ ^ | X) = cos 9 \ X) — i sin 9\Y). 
The (T^ — rotation is performed at the effective spin 
that connects the pair of X particles. A relative phase 
between the two states | Ol) and can be generated 
by changing the relative self-energy of the correspond- 
ing excitations. Hence, a general one-qubit gate can be 
obtained in the usual dynamical way. 

The realisation of a two qubit gate requires braid- 
ing of the anyonic particles. To transport the anyons 
around the two dimensional lattice one can create a trap- 
ping potential well with depth V for each anyon. This 
can be achieved by reducing the coupling Jcs by V at 
the corresponding hexagon. A simple implementation 
is by increasing the couplings of the two z-links that 
surround a particular hexagon occupied by the anyon 
to the value J^. Thus, the effective coupling becomes 
= Jjj2/(4(J^ -1-4)2 J^). If Ji = 3J^, then the eigen- 
value of the excitation corresponding to that hexagon 
is reduced twelvefold generating a strong and efficient 
trapping well. The employed interaction, a^a^ , does not 
change the nature of the particles, i.e. it does not cor- 
respond to a a^, CT^ or rotations. By moving the 
trapping potential in an adiabatic fashion it is possible 
to move each excitation independently. Hence, braiding 
can be performed that can actually reveal the anyonic 
behaviour of the particles or can be employed to process 
information. 

Now we have all the necessary ingredients to build a 
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quantum computer. The register includes pairs of parti- 
cles so one can arrange them as in Figure |21 Each qubit 
is allocated at a specific region on the lattice where ar- 
bitrary superpositions of X- and y-particles can be pro- 
duced to encode an arbitrary qubit state. To perform a 
two-qubit gate one needs to transport the X-particle of 
the control qubit around both, X- and y-particles of the 
target qubit, as shown in Figure El As a result the state 
I XY) acquires a minus sign while all the other states re- 
main te same. In the logical space this corresponds to 
a controlled phase gate, where a minus sign is generated 
only for the logical state | 01). This procedure can be per- 
formed between any two arbitrary qubits, which together 
with the general one-qubit gate result into universality. 



qubit 1 



qubit 2 



qubit 3 



FIG. 2: A possible arrangement of anyons comprising a quan- 
tum computer. The red particles are the fermions, X, while 
the blue ones can be either Y or Z particles. The depicted 
evolution results in a controlled phase gate. 
The value of the controlled phase gate is fixed by the 
statistics of the permuted particles and is resilient to any 
deformations of the spanned loop as long as it circulates 
the target particles once. Moreover, moving the parti- 
cles around with the trapping potentials is immune to 
control errors that would occur if the transportation was 
performed by dynamical means. For the latter case one 
has to rotate the eS'ective spins successively along the 
desired trajectory, with errors that would accumulate at 
each translation step. By moving the potential well adi- 
abatically the trapped particle follows the lowest energy 
configuration without the need for any dynamical control 
of the effective spins and subsequent introduction of con- 
trol errors. This resilience which is of geometrical origin 
gives a great advantage in transporting and braiding the 
anyons fault tolerantly. 

It has been shown [To] that Hamiltonian ^ can be gen- 
erated with a Bose-Einstein condensate of atoms super- 
posed with a two dimensional optical lattice configura- 
tion. In particular, consider two bosonic species labelled 
by (7 = a, 6, e.g. given by two hyperfine levels of an atom. 
They can be trapped by employing two in-phase optical 
lattices for each direction. The geometry of the hexago- 
nal lattice can be produced with three pairs of properly 
aligned counterpropagating laser fields ^^l- The tun- 
nelling of atoms between neighboring sites is described by 
— X]!o-(*cra/crO(z+i)cr+H.c.). Whcu two or more atoms are 
present in the same site, they experience collisions given 
by i U a\^a\^, ai^i aia ■ We now consider the limit 



t ^ U where the system is in the Mott insulator regime 
with one atom per lattice site In this regime, we can 
employ the pseudospin basis of |t) = jn^* = 1, nj* = 0) and 
J.) = = 0,rtf = 1), for lattice site I, and the effective 
evolution can be expressed in terms of the correspond- 
ing Pauli operators. Consider initially high enough laser 
intensities such that each site corresponds to a free spin 
1 /2 particle. By applying Raman transitions between the 
states I a) and | b) it is possible to perform an arbitrary lo- 
cal spin rotation. It is easily verified that when one of the 
tunnelling coupling is activated, say ta, the Ising inter- 
action is realized between neighboring sites 0, 0, 0| , 
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where < l,m > are the nearest neighbors connected by 
the ta coupling. This procedure can implement the in- 
teractions of the 2-links of Hamiltonian QJ, where the 
Jz coupling varies as a function of ta- To generate the 
interactions in the x- and y-links one needs to create 
an anisotropy between the x and y spin directions. For 
that we activate a tunnelling by means of Raman cou- 
pling consisting of two standing lasers Li and 
L2, giving finally the interaction terms J2<i m> '^f'^m '^^ 
S<z m> '^I'^m depending on the phase difference between 
the laser radiations, Li and Appropriate combina- 
tions of these interactions gives finally all the necessary 
terms for realizing Kitaev's Hamiltonian. 

Upon this basis we can realize a quantum computa- 
tional scheme with abelian anyons. As we have seen, the 
generation of particle excitations involves the rotation of 
effective spins, or equivalently, rotations of the spins at 
the z-links of the original lattice. In this way initializa- 
tion of the qubit space is straightforwardly implemented 
by local Raman transitions. To perform a one-qubit rota- 
tion an appropriate detuned laser radiation can be used 
to induce a local rotation. To perform the trapping 
of the excitations we choose to use localized a^cr^ inter- 
actions that would lower the Jeff coupling, thereby re- 
ducing the overall energy of the excitations at a specific 
plaquette. This can be implemented simply by focusing 
a laser field in-between the spins of the 2:-link, that will 
reduce the potential barrier thus increasing the ta tun- 
nelling (see Eqn. Note, that this laser does not 
need to have a cross-section smaller than the lattice pe- 
riod, but merely to have a strong effect on the tunnelling 
between two certain sites and weaker between the rest 
of the sites, thus trapping the anyonic excitation in that 
neighborhood. By moving the focus of this laser around 
the lattice, the potential minimum and its trapped parti- 
cle will also move accordingly. To suppress the generation 
of undesired relative phases between the excitations due 
to fiuctuations of the intensity of the trapping laser one 
can use the same laser source for generating the trapping 
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potentials for all the particles. The main source of de- 
coherence, beyond the leakage of coherence outside the 
pseudospin space due to heating of the atoms, is the gen- 
eration of virtual pairs of particles due to thermal fluctu- 
ations. These excitations can actually move the particles 
out of their trapping potentials, which will register as 
decoherence of the encoded information. To avoid such 
problems, we would like to have the temperature small 
enough compared to the depth of the trapping well, V, 
thus exponentially suppressing the generation of virtual 
pairs. 




FIG. 3: (a) The numerically obtained energy gap, AE, be- 
tween the ground and the first excited state with Jz = 1. (b) 
Loglog plot of AE as a function of the coupling J = Jx ~ Jy 
There is very good agreement between the solid (blue) line, 
predicted by perturbation theory, and the dashed (red) line, 
given by numerical diagonalization. (c) AE, for various states 
above the ground state. As expected from the theory, AE al- 
most doubles for a certain set of states that correspond to X 
excitations compared to Y or Z excitation. 

The measurement procedure required for the presented 
scheme has to be able to distinguish between the vacuum 
state and a particular excitation. As we have seen, exci- 
tations are localized objects that live within the bound- 
aries of the trapping potentials. A generic measurement 
should therefore be able to address the density profile 
of the particular atomic states as a function of position 
and determine their correlations. To increase the distin- 
guishability between their correlations one may alterna- 
tively use local time modulations of the optical lattice to 
induce Bragg scattering resonant between the particular 
excitation and the first excited vibrational state of the 
atoms. By direct population measurement or by mea- 
surement of the atomic correlations one can, in principle, 
deduce the state of the system [H E [H . 

Finally, we would like to see to which extend pertur- 
bation theory describes the initial model by employing 
numerical diagonalization of Hamiltonian Towards 
that we consider a periodic honeycomb lattice consist- 
ing of 16 spins. Figure EJa) presents the energy gap. 



AE, between the ground and the first excited states as 
a function of Jx and Jy where we have set = I. It 
is clear that the model has a gap for |Ja;| + I'^yl < \Jz\ 
as expected from the analytical treatment 0. Further- 
more, it is possible to compare the energy gap of the 
perturbative theory with the one obtained numerically 
for Jx,Jy Jz- Figure I^Jb) presents a comparison be- 
tween perturbation (solid, blue line) and the numerical 
diagonalization (dashed, red line) of AE for various val- 
ues of the coupling J — Jx = Jy- We observe that the 
two curves are in excellent agreement indicating the va- 
lidity of perturbation theory for a wide range of coupling 
values J. As a last test we consider a part of the spec- 
trum above the ground state. We saw previously that the 
energy gap of the fermionic excitation, X, is double than 
the energy gap of F or Z excitation. Indeed, Figure|2Ic) 
gives AE for the first few excitations where the doubling 
of AE is clearly observed. Unfortunately, the system of 
16 spins that we are able to diagonalize numerically is too 
small to represent faithfully the properties of the spec- 
trum expected for a large system. A complete theoretical 
study that presents the properties of the excitation be- 
yond perturbation theory will be presented elsewhere. 
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